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Description 

[0001 ] The invention relates to a method for determin- 
ing a contour in a space with a density variation, com- 
prising: 5 

defining a plurality of initial vertices being located 
on a seed contour, said vertices being connected 
by edges to neighbouring vertices; 
defining ah energy function comprising an internal 10 
portion which is a function of the curvature at each 
of the vertices and an external portion determined 
by the density variation at the location of each of 
said vertices; 

determining a final contour by continuously varying 15 
the position of the vertices till the energy function 
reaches a minimum. 

Such a method can, for example, be used in a medical 
environment for (semi-)automatically defining the out- 20 
line of a region of interest or of an object, such as an 
organ or a tumour, in a 2 or 3 dimensional image. Images 
of clinical importance can be obtained by a variety of 
modalities for example, a conventional X-ray technique, 
a CT-scanner, magnetic resonance imaging, single pho- 25 
ton emission computerised tomography, positron emis- 
sion tomography or by means of ultrasound echoscopy. 
Other areas in which deformable contours can be used 
are, for example, computer graphics and animation. The 
defined outline of the object can, for example, be used 30 
as a basis for providing quantitative information to a phy- 
sician, for surface extraction for visualisation or for vol- 
ume definition. 

[0002] The term energy function is not related to phys- 
ical energy, but is used because the procedure of con- 35 
tour definition has some resemblance to a string of in- 
terconnected masses moving in a field of varying poten- 
tial energy and trying to establish a stable situation of 
minimum energy. In this similarity the external portion of 
the energy function corresponds to the density or poten- *o 
tial energy and the internal portion to the mutual inter- 
action of the elements of the string. 
[0003] Such a method for defining is known from an 
article by J.V. Miller et a/., entitled "Geometrically de- 
formed models: A method for extracting closed geomet- 45 
ric models from volume data", published in Computer 
Graphics, Vol. 25, No. 4, (1991), pages 217-226. In that 
article a contour is described as a set of vertices con- 
nected by edges. The energy function comprises a to- 
pology preserving energy term, dependent on an esti- 50 
mation of local curvature and the distance between a 
vertex and its neighbours, an image event energy term, 
derived from the density (or pixel values), and a locally 
defined deformation potential driving the vertices out- 
ward or inward. The energy function is evaluated for the 55 
vertex positions, not for the trajectory of the connecting 
edges. This makes the contour discrete, whereby the 
resolution is determined by the length of the edges. 



Processing of a seed contour, entered by an operator 
as a number of vertex points, to a final contour occurs 
in an iterative procedure. During each step the vertices 
are moved in the direction of steepest descent along the 
surface provided by the energy function. Movement of 
a vertex stops when no energy reduction occurs. 
[0004] In the known method, the value of the energy 
function can decrease if a vertex is displaced along an 
edge and, therefore, the different vertices tend to cluster 
in corners of the contour. Another disadvantage is that 
the contour may collapse or expand within a region in 
which no density variation is present. To avoid such be- 
haviour the known method needs the topology preserv- 
ing term in the energy function. This topology maintain- 
ing term provides constraints or interactions between 
vertices which are artificial and not related to the density 
variation in the space. Accordingly, deviations between 
the contour obtained and the actual shape of the object 
may occur due to these constraints. 
[0005] It is, inter alia, an object of the invention to pro- 
vide a method for defining a contour in which method 
there is no interaction between vertices. 
[0006] This object is achieved by the method accord- 
ing to the invention as defined in Claim 1 . Another object 
of the invention to provide an arrangement to carry out 
the method of the invention is achieved by the arrange- 
ment as defined in Claim 9. The energy contribution of 
each vertex to the contour as a whole does not depend 
on the distance to the neighbouring vertices and there 
is no direct interaction between the vertices. As the en- 
ergy variation at a vertex depends only on the direction 
"perpendicular" to the local path of the contour, during 
each iteration vertices will move only in this perpendic- 
ular direction and clustering of vertices will be avoided. 
[0007] The method according to the invention can be 
further characterised in that the internal energy portion 
is derived from the difference in angles of the edges at 
a series of adjacent vertices. In order to avoid shrinking 
purely due to the internal portion of the energy function, 
the internal portion is made dependent on the shape of 
the contour. In this embodiment the internal portion will 
only contribute if the "curvature" of the contour, i.e. the 
numerical value of the angle between edges at adjacent 
vertices, changes. When the internal portion of the en- 
ergy function contributes in this way, in a region with no 
density variation a circular portion of the contour will not 
shrink, while a zigzag portion of the contour will straiten. 
Both kinds of behaviour represent a desired feature of 
the contour. 

[0008] An embodiment of this method is character- 
ised in that said internal energy portion is obtained by 
the convolution of said difference in angles with a sym- 
metric discrete filter having a zero frequency component 
equal to zero. Such a filter imposes the desired behav- 
iour. In a two-dimensional contour the most simple form 
of such filter is (-2, 1, -2), i.e. the internal portion at a 
vertex is proportional to the angle at the vertex itself mi- 
nus the average of the angles at both adjacent vertices. 
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[0009] An embodiment of the method according to the 
invention is characterised in that the seed contour is 
open. In a two-dimensional situation this corresponds to 
a curve in the space having two ends. Such an end has 
only one edge and no energy contribution from a term 5 
that corresponds to the angle of two edges. In case of 
a three-dimensional contour it means that a face be- 
tween a number of vertices is, at least at one side, not 
joined by another face. In both cases vertices can move 
away from the contour, thereby extending the contour 
and eventually closing it 

[0010] A preferred embodiment of the method accord- 
ing to the invention is characterised in that if during min- 
imization of the energy function an edge length exceeds 
a first predetermined threshold, the edge is replaced by 
an additional vertex and connecting edges. Although, in 
the present method, the movement of a vertex along an 
edge will not change the energy contribution, vertices 
may move away from each other due to energy contri- 
bution of the density variation. Too long edges reduce 
the resolution of the contour. To maintain the resolution 
at regular intervals, vertices are added in long edges. 
[001 1 ] The opposite may also occur, if edges become 
too short a computational load in performing the method 
is present without any benefit in resolution. In that case 
the short edge is taken out. 

[0012] An embodiment of the method is characterised 
in that a user defined energy distribution is added to the 
external portion- of the energy function. Such a user 
function forces the contour to respect user defined re- 
strictions, for example to keep out of areas in the space 
where the user knows no interesting features are 
present. 

[0013] If contours have to be determined in a number 
of similar images, for example images representing ad- 
jacent slices in a body obtained by CT or MRI, it is pos- 
sible to use the results in one slice as a starting point in 
the next. Such an embodiment of the method is charac- 
terised in that firstly a contour is determined in a first 
space and that, secondly, a contour is determined in a 
second space having a density variation having similar 
features to the density variation in the first space where- 
by the seed contour in the second space is defined by 
transfer of the final contour obtained in the first space. 
[0014] The invention also relates to an arrangement 
for determining a contour a space in which a given pa- 
rameter has a value varying throughout said space, the 
arrangement comprising 

a memory for storing parameter values in said 
space in the form of a matrix of pixel values; 
a display having a two-dimensional matrix of pixels; 
pointer means for indicating initial vertex positions 
for the definition of a seed contour; 
- computing means, having access to said memory 
and said vertex positions, for evaluating an energy 
function comprising an internal portion which is a 
function of the curvature at each of the vertices and 



an external portion determined by the density vari- 
ation at the location of each of said vertices; 

the computing means being arranged for 

determining a final contour by continuously varying 
the position of the vertices till the energy function 
reaches a minimum, while 
for each of the vertices the internal portion of the 
energy function is dependent on the angle between 
the edges at the respective vertex and in that the 
external portion of the energy function is related to 
the density variation in a direction having equal an- 
gles with each of the edges at the respective vertex. 

Such an arrangement or workstation is used to perform 

the method according to the invention. 

[001 5] These, and other more detailed aspects of the 

invention will now be elucidated by way of example with 

reference to the accompanying drawings. 

[0016] The drawings show in 

Figure 1 the basic structure of contour and its ele- 
ments as used in the method according to the 
present invention; 

Figure 2 an illustration of the internal forces and 
their effect on the contour; 
Figure 3 an illustration of the definition of local cur- 
vature at a vertex in the method according to the 
invention; 

Figure 4 an illustration of the radial and tangential 
directions in a vertex of a contour; 
Figures 5a, 5b and 5c, the effects of uncompensat- 
ed internal forces in typical contour shapes and in 
absence of external forces; 
Figures 6a, 6b and 6c, the internal forces in typical 
contour when derived by means of a filter taking into 
account the variation from vertex to vertex; 
Figure 7 the locally radial components f im „• of the 
external forces f jm v/ according to the invention that 
provide a resulting driving force on the vertices of 
the contour; 

Figures 8a and 8b the replacement of two nearby 
vertices by a single vertex and the insertion of an 
additional vertex in between adjacent vertices that 
are far away from each other; 
Figure 9 diagrammatically a workstation for per- 
forming the invented method. 

[0017] For reasons of convenience, in the following 
description the effects of the minimisation of the energy 
function are described in the term of forces acting on the 
vertices. In the similarity between a string of particles or 
masses moving in a potential energy field the two ap- 
proaches are equivalent in that a particle in an energy 
field can be thought of as being subject to a force that 
is related to the gradient of the energy: F = -VE. The two 
ways of description are completely equivalent. A de- 
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scription in terms of forces however,. seems to be more 
convenient for the present method, because the forces 
relate directly to the acceleration and displacement of 
the vertices of the contour. 

[0018] In the description the method is illustrated by 5 
a two-dimensional version. Extension to three-dimen- 
sional spaces is straightforward, in as much as the struc- 
ture will be edge-connected vertices and evaluation will 
be performed locally on vertex positions. 
[0019] Starting from an initial shape, which can be *0 
created with a minimum of user interaction, the dynamic 
contour actively modifies its shape, thus approximating 
some desired contour. The driving force behind the 
shape deformation is calculated from internal forces, de- 
rived from the shape of the contour itself, and an exter- 15 
nal force field, derived from some image feature energy 
distribution. The internal forces will try to minimize local 
contour curvature, while the external forces will try to 
make the contour follow a valley or a ridge through the 
'landscape' formed by the image feature or density var- 20 
iation. By applying both internal and external forces with 
user definable weight factors, the user can determine to 
either follow the image feature landscape in a very glo- 
bal way, or very precisely, or anything in-between. 
[0020] The deformation process is performed in a 25 
number of discrete steps, after each of which the situa- 
tion with regard to position, velocity and acceleration is 
evaluated for each of the vertices. In this evaluation in- 
ternal and external forces on a vertex are calculated 
from the position of the vertex and its neighbours. These 30 
forces result in an acceleration, which changes the ve- 
locity of the vertex. This velocity determines the dis- 
placement of the vertex during the next deformation 
step. After a number of deformation steps a stable end 
situation will be reached in which there is an equilibrium, 35 
which means that velocity and acceleration are zero for 
each vertex. Described in terms of energies, this situa- 
tion represents a local minimum of the energy function. 
[0021] During deformation there are two undesirable 
effects which may occur: shrinking of closed contours *o 
owing to internal forces, and clustering or gathering of 
vertices in corners of the contour, owing to external forc- 
es. In the present invention the shrinking problem is 
overcome by a proper definition of local curvature at the 
vertices, combined with a definition of internal forces 45 
which is derived from local curvature in such a way, that 
these forces are not only zero for parts of the contour 
where curvature is zero, but also for parts where the cur- 
vature is constant. The solution to the second problem 
is found by allowing only locally radial vertex displace- 50 
ments. 

[0022] Figure 1 presents the basic structure of the 
contour 1 0 as used in the invented method. The contour 
is an object consisting of vertices 11, 12, 13, ... which 
are connected by straight line segments or edges 14, 55 
15, 16, .... The position of a vertex V, is represented by 
a vector p, and the edge between V t and V i+1 by a vector 
d/, the coordinate system is assumed to be Cartesian. 



Deformation is caused by a combination of forces which 
act on the vertices; the resulting acceleration in vertex 
Vj is denoted by a vector a,-. Another property of a vertex, 
not shown in the Figure, but important for the dynamic 
behaviour of the contour, is its velocity, denoted by v, for 
vertex V,. 

[0023] As mentioned before, the length ||d,|| of an edge 
segment represents the local resolution of the contour: 
if it is large the contour will not be able to follow varia- 
tions of small scale in the image feature energy distri- 
bution. The length ||d J may change during every defor- 
mation step, thus causing local variation in the resolu- 
tion of the contour. To keep this variation limited, the 
edge lengths are evaluated at regular intervals; where 
necessary, vertices are removed or inserted, thus keep- 
ing the resolution of the contour close to a user-specified 
scale. 

[0024] The energy function comprises an internal por- 
tion or, in other words, the vertices are subject to internal 
forces. The internal forces defined in the present meth- 
od are related to the local curvature of the contour. The 
main objective of introducing internal forces or energy 
functions is illustrated in Figure 2: minimizing local cur- 
vature, forming a counterbalance to external forces that 
try to shape the contour according to all the variations 
of the image feature landscape, in order to obtain a good 
approximation of a smooth curve 20. 
[0025] Firstly, the concept of local curvature in a dis- 
crete contour is to be defined, which is not a trivial mat- 
ter. Strictly speaking, local curvature is zero on the 
straight edge segments in-between the vertices, while 
it is not defined at the exact position of a vertex (first 
order discontinuity), which however happens to be ex- 
actly the position where it needs to be defined. A satis- 
factory solution is defining local curvature at the position 
of a vertex to be the difference between the directions 
of the two edge segments that join at that location. The 
edge segment leaving from vertex V, is represented by 
a. vector d,; its direction is described by the unit vector 
d/. According to the definition above, the local curvature 
c,at V, is described by (see Figure 3): C/= d r d M . Defined 
in this way, local curvature has length (strength) as well 
as direction and provides a usable and unique measure 
for the angle between two joining edge segments. More- 
over, the length of the curvature vector depends only on 
this angle, and is not influenced by the lengths of the 
two joining edge segments. 

[0026] Also locally radial and tangential directions at 
the position of a vertex are defined. For this use is made 
of the unit vectors d, representing the directions of the 
edge segments d,-.Jn the present method the locally tan- 
gential unit vector t, is defined as the normalized sum of 
the unit vefctojs of two/joining edge segments (see Fig- 
ure 4): t, = d r d M / ||d,+d M ||. Thermit vector 17 in the local 
radial direction is derjyed frgm t y by a rotation over nJ2 
radians. The vectors t, and 17 now represent a local co- 
ordinate system at the position of vertex V h which is use- 
ful for the calculation of internal as well as external forc- 
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es. 

[0027] Both closed and open contours are allowed. If 
the number of vertices equals n and the contour is 
closed, the first and the last vertices are connected, so 
V 0 has two neighbours: V n _ 1 and V v If, however, the 5 
contour is open, V 0 and V n _ A are not connected and have 
both but one neighbour: V Q only connects to and v' n _ 1 
only connects to V n _ 2 . This situation requires special 
measures for the, calculation of the local tangential and 
radial directions t,- and r, as well as the curvature vector 10 
Cj. At the position of the open ends the local tangential 
direction is taken to be equal to the direction of the first 
c.q. the last contour segment. The length of the curva- 
ture vector is set to zero for both end positions. 
[0028] The local curvature vector c,, when described 15 
in terms of the local r,f-coordinate system is always di- 
rected along the^ local r-axis, pointing either in exactly 
the direction of r h or in the opposite direction. In other 
words, Cj is a vector along the local r-axjs and its length 
can be described by the dot product (c,T;). According to 20 
this definition, the length ||cj| of the curva^tu/e vector can 
be positive as well as negative: c,- = (C/T,)^. 
[0029] Local curvature is defined now as a one-di- 
mensional variable in locally radial direction. The inter- 
nal forces which should act on the vertices of the contour 25 
are to be defined arid restrictions put on the deformation 
process. In order to get a clear understanding of the con- 
tribution that internal forces will make to the deformation 
of the contour, a situation is considered in which external 
forces are completely absent. 30 
[0030] As can be seen from the Figures 2 and 3, the 
desired internal forces (Figure 2) and the curvature vec- 
tors (Figure 3) are very closely related; both are vectors 
in locally radial direction with the same orientation. It 
would however not be wise to define the internal forces 35 
as being proportional to the local curvature vectors. The 
reason for this is illustrated in the left part of Figures 5a, 
5b and 5c. Any simple closed shape, such as shown in 
Figure 5a in Cartesian representation, would in the ab- 
sence of external forces be deformed into the shape 21 *o 
having the minimum overall local curvature, i.e. a circle 
(or rather, because of the discretisation, a symmetric 
polygon). Then however, the deformation process 
would not stop, but continue to move the vertices in the 
direction of the centre of the contour. During this last 45 
phase the contour is shrinking while local curvature will 
not change, which means that the internal forces are not 
doing what they were intended to do: reduce local cur- 
vature. The contour will completely implode into a single 
point. 50 
[0031] By introducing a second constraint on the po- 
sition of a vertex with respect to its neighbours, the 
shrinking process can be brought to a halt. This second 
constraint acts like an elastic force that keeps the dis- 
tance between neighbouring vertices between certain 55 
limits. Vertices which are moving closer together will ex- 
perience an increasingly repelling force that at some 
point will stop the shrinking process. The point where 



the shrinking process comes to a stop depends on the 
weights that are assigned to the two internal constraint 
forces. In other words: it is required to optimize the bal- 
ance between the internal elastic force and the internal 
curvature minimizing force. 

[0032] In order to avoid the introduction of the second 
- rather artificial - constraint and the problems entailed 
by it, in the method according to the present invention it 
is required that local curvature should be reduced with- 
out affecting parts of the contour with constant curva- 
ture. 

[0033] Figures 5a, 5b and 5c illustrates on the left in 
Cartesian representations, besides the already dis- 
cussed symmetric polygon 2 1 in Figure 5a, some typical 
shapes which may occur in a discrete contour. Figure 
5b shows a portion 22 of a contour having a turn in di- 
rection over k radians and Figure 5c a portion 23 of a 
contour with alternating curvature direction. This further 
illustrates the inadequacy of using an internal force pro- 
portional to the curvature vector. 
[0034] The shape 22 in Figure 5b also shows a prob- 
lem. It forms an extended part of a contour. If a force 
proportional to local curvature would be applied to the 
vertices, the extended part would become shorter as 
shown in the Figure. However, as with the shrinking of 
the shape in Figure 5a, this would not actually reduce 
local curvature, but only displace the curved region. 
[0035] The shape 23 shown in Figure 5c does not 
pose a problem; it is included because it represents a 
typical situation in which local curvature reduction may 
be required. Any solution to the shrinking problem 
should still perform well on this shape. 
[0036] On the right side of Figures 5a, 5b and 5c the 
same three shapes are shown, but now in local ^-coor- 
dinates, labelled 21\ 22' and 23', respectively. Exami- 
nation of the curvature vectors in this coordinate system, 
combined with the intention to derive from these vectors 
the internal forces that would meet the above stated cri- 
teria, provides a solution. First: internal forces f in / which 
act on the vertices V, should have the same (radial) di- 
rection as the curvature vectors. This means that inter- 
nal forces can be derived from the curvature vectors by 
modifying only their lengths. Second: in order to reduce 
local curvature without affecting areas of constant cur- 
vature, the lengths of the internal force vectors should 
be zero for parts of the contour with constant curvature. 
Both conditions can be met if the sequence of dot prod- 
ucts (CjT/) along the contour is considered as a discrete 
scalar function, and if the convolution of this function 
with a discrete filter /c, is used as the representation of 
the A sequence of internal force vector lengths ||f ir J| = 

[0037] The first condition cafl be met by using r,- as 
the direction of f in / : f ln/ = ||f jr J| - r,-. The second condition 
is met by choosing appropriate filter coefficients for k f in 
which a convolution of a sequence of constant values 
with kj will be zero. This is achieved by a filter /c, which 
is a symmetric discrete filter with a zero frequency com- 
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ponent equal to zero. A wide class of filters complies 
with this condition. In its simplest form it is a filter con- 
sisting of three coefficients with the values k f = 1 , 
- 1 /£). Other possible filters provide an interesting way to 
optimize the behaviour of the contour for specific appli- 5 
cations. In particular, in may be advantageous to make 
kj adaptive. 

[0038] The results of this filter are illustrated in Fig- 
ures 6a, 6b and 6c, which shows on the left side again 
the three typical shapes 2 V, 22' and 23' in local r,f-co- 10 
ordinates, but now with the internal force vectors which 
were derived from the curvature vectors. The right side 
of the Figure shows the shapes 21, 22 and 23 and in- 
ternal forces in Cartesian coordinates. The constant cur- 
vature of shape 21 in Figure 6a of the Figure resulted in 15 
internal forces that are zero everywhere, showing that 
the shrinking problem has been solved. The n radians 
turn of the shape 22 of Figure 6b is deformed in a way 
that appears more natural. Not only does the extended 
contour part become shorter, but it also becomes wider, 20 
which is what would be expected from curvature reduc- 
tion. The alternating curvature 23, shown in Figure 6c, 
will still be reduced effectively, as the effect of convolu- 
tion with kj on a purely alternating signal comes down to 
multiplication with a constant. In consequence, -the in- 25 
ternal forces f in ,, as defined hereinbefore, produced the 
desired deformation effects to the contour in the ab- 
sence of any external forces. 

[0039] To provide a driving force for the deformation 
of the contour, an external potential energy distribution 30 
is assumed, which represents the energy or strength of 
some kind of image feature or combination of image fea- 
tures. The selection of appropriate image features for 
specific applications is an important subject but not part 
of the invention. In the following the behaviour of the 35 
contour using simple image features like the pixel or 
voxel grayvalue itself and the length of the grayvalue 
gradient is demonstrated. If, for instance, a contour 
should follow a maximum gradient path through the im- 
age, the gradient length could be used as an image fea- 40 
ture and an energy distribution defined that is high for 
large values of this feature. A maximum gradient path 
is then represented by a ridge in this energy distribution. 
The implementation of the deformation process is such 
that it will attempt to pull the vertices into local minima 45 
of the energy distribution, which is a natural behaviour. 
For the contour as a whole this means that it will end up 
following a path of low energy or a valley through the 
feature energy landscape. If the contour is to follow a 
ridge instead of a valley, as in this example, this can be so 
done by inverting the energy distribution. The resulting 
distribution of potential energy is called E im . The force 
field that will pull an object in the direction of lower en- 
ergy, can be described by the following simple relation: 

[0040] If this force is applied to the vertices of a con- 
tour in a situation in which there are no internal forces, 
then the result will be that in the end the contour con- 



nects local energy minima, following a valley through the 
external energy distribution. However, in case of apply- 
ing the above described force field, a problem is that the 
force f jm , which acts on the vertex V h will in general have 
a component along the path of the contour (locally tan- 
gential component). This tangential component may be 
substantial, or even dominant. If there were no restric- 
tions on the curvature of the contour (e.g. in the absence 
of internal forces), then in the final situation the contour 
would actually pass through the local minima of the ex- 
ternal energy distribution, in which case the force field 
would locally be directed along the path of the contour! 
The result of this locally tangential component is that 
vertices will move along the contour and form clusters 
in local minima of the external energy distribution, which 
is obviously a very undesirable effect. 
[0041] Introducing an elastic force as a second inter- 
nal force, in order to keep the distance between neigh- 
bouring vertices within limits, assuming that such a force 
could be an effective answer to the clustering problem. 
In this situation the strength of the elastic force would 
have to be tuned locally to the strength of the external 
force f im Vi at the location of vertex V f . This external force 
varies, not only over the area of the image matrix, but 
also from image to image and it depends on many pa- 
rameters related to both data acquisition and processing 
of the acquired data, like for instance the feature extrac- 
tion algorithm that was used. Careful tuning of the inter- 
nal elastic force is important: If the elastic force is too 
weak it may lead to shrinking of the contour owing to 
internal curvature forces, or it may lead to clustering of 
the vertices owing to external forces. If the elastic force 
is too strong it may obstruct the deformation process of 
the contour, because during deformation the vertices 
must have the freedom to move and change the dis- 
tance to their neighbours. It follows that correct tuning 
is very important, but in general not possible. The inter- 
nal elastic force cannot be tuned locally to the external 
forces f im Vi and at the same time be tuned correctly to 
the internal forces f jn 

[0042] According to the present invention, this prob- 
lem is solved in line with the solution to the shrinking 
problem, discussed hereinbefore. Vertex displacement 
along the path of the contour does not make any contri- 
bution to the deformation. Only the locally radial com- 
ponent of f im Vi is used to drive the vertices of the con- 
tour. If the locally radial component f im Vj is denoted by 
f im f7 -, A then its length js given by the dot product of f im Vi 
and r,: f jm „• = (fj mt w" r /) r /- Tnis is illustrated in Figure 7, 
which shows the locally radial forces f im r/ and f jn , that 
provide a resulting driving force on the vertices of the 
contour which is purely devoted to deformation of the 
contour, without moving vertices along the path of the 
contour. 

[0043] For a user of dynamic contour modelling tools 
to have control over the contouring process is an impor- 
tant aspect. Normally, real world images lend them- 
selves rarely for fully automatic contouring to give sat- 
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isfactory quality. Some operator guidance may be desir- 
able for object definition purposes. Also, the state of 
minimum energy reached by the contour may often be 
one out of a set of possible-local minima, and the oper- 
ator should then have the possibility to push the contour 5 
out of a local minimum into another. Therefore, the 
present method, as described hereinbefore, may be 
supplemented with a possibility for the user to interac- 
tively create and modify a user defined external energy 
distribution E user , to be added to the calculated image 
feature energy distribution E jm to form a combined ex- 
ternal energy distribution E ex . The user defined distribu- 
tion enables the operator to create additional ridges and 
valleys through the energy landscape to locally force the 
contour to follow a particular path. Once these two ex- 
ternal energy distributions are combined, the situation 
is completely equivalent to the one above, and calcula- 
tion of the resulting external force f ex r/l action a ver- 
tex V,of the contour becomes: f ex r/ = (f e x,w r /) r /, ' n which 

fex=-V(£im + E user ). 

[0044] With these definitions of the internal and exter- 
nal forces used in the method the dynamics of the de- 
formation process according to the method will be de- 
scribed. 

[0045] The distance between a vertex and its neigh- 
bours determines the resolution of the contour, details 
of the external image feature energy distribution, which 
are small enough, may pass through the spaces be- 
tween the vertices, without having a significant influence 
on the final shape of the contour. Because of the defor- 
mation of the contour, the vertex to vertex distance will 
change constantly. This may result in local variation as 
well as in global changes in the resolution of the contour. 
Both are undesirable, but slight local variation in the res- 
olution is unavoidable, because a deformation process 
can only be possible if the vertices have the freedom to 
move with respect to their neighbours. For this reason, 
the variation is kept between certain limits, by periodi- 
cally resampling the contour along its path. In local r,t- 
coordinates this translates into a resampling along the 
f-axis. 

[0046] The total force f, acting on a vertex is a weight- 
ed combination of external and internal forces: f, = 
w ex f ex r/ + w in f jn/ . The weighting factors w ex and w in may 
have default values for each application but allow mod- 
ification by the operator. Upon emphasizing on external 
forces one will make the contour follow the extracted im- 
age features more precisely, while putting emphasis on 
internal forces will smooth out the path of the contour. 
As a result of the forces that act on a vertex V h this vertex so 
will start to move and change its position p,. This position 
vector, together with the vertex velocity and acceleration 
vectors v, and a,-, describe the dynamic state of a vertex. 
A vertex will not stop moving until both ||vj| = 0 and ||aj| 
= 0. The deformation process of a contour is not com- 55 
pleted before this condition ||vj| = ||aj| = 0 is met by all its 
vertices. It may in principle take very long before the 
contour comes to a rest, or the contour may even remain 



oscillating between two states, which both represent a 
local energy minimum. To avoid such behaviour a third 
component is added to the force that is applied to a ver- 
tex V h viz. a damping force, proportional to the vertex 
velocity v,: f, = w ex f QXJi + w in f jn / + w darnp v,,The weighting 
factor w damp is negative and determines the amount of 
damping. Even a small damping factor suffices to en- 
sure stability of the deformation process. 
[0047] In an embodiment the actual deformation proc- 
ess is implemented as a so-called numerical time inte- 
gration process, in which the complete state of the con- 
tour is calculated at a sequence of discrete positions in 
time. Let the argument t describe the state of the contour 
at a certain point t in time, and r+Af the situation at a 
time At later, then the deformation process during the 
incremental time At can be described by the following 
difference scheme: 

P/ (f+A0 = P / (r) + v / (f)Af; 
v / (r+A0 = v / (f) + a / (0Ar; 

a^f+Af) = f,(f + At)lm r 

The value of f,{f + At) is calculated as indicated herein- 
before. The value m, would, in a physical context, rep- 
resent the 'mass' of vertex V h In a prototype implemen- 
tation of the method no use is made of the possibility to 
assign different mass values to the vertices of the con- 
tour. By assuming the same mass for all vertices, m t re- 
duces to a constant scaling factor. In the implementation 
the procedure was further simplified by giving the time 
interval At a value of one. 

[0048] An option which is available for open contours 
is a growing process that adds vertices to one or both 
of the open ends of the contour. By integrating this grow- 
ing with the deformation of the contour an efficient track- 
ing method is obtained. Vertices that are added are po- 
sitioned by extrapolating in the direction of the last (or 
first) edge segment, which is at the same time continu- 
ously being repositioned by the deformation process. An 
automatic closing option can be added as well, such that 
both ends of an open contour are automatically connect- 
ed when the distance between them becomes smaller 
than a certain threshold. 

[0049] In a prototype that was implemented, the user 
can interactively control a parameter / des , the desired 
length for the edge segments of the contour. In this way, 
the resolution of the contour can be chosen. A default 
value for / des can be calculated by normalizing / des by 
the inner scale of the image. This means that / des in- 
creases when the image is considered at a coarser 
scale. From / des two other values / min and / max are de- 
rived, representing the minimum and maximum dis- 
tance which is allowed between neighbouring vertices. 



15 



20 



25 



30 



35 



40 



7 



13 



EP 0 633 548 B1 



14 



[0050] In the prototype, resampling is implemented as 
a two pass process: The first pass checks along the en- 
tire contour if any segment length has become shorter 
than the minimum length / min . If such an edge segment 
is found, this edge segment is removed from the contour 5 
by replacing the two vertices on both ends of this seg- 
ment by one single vertex at a position exactly in-be- 
tween the replaced vertices. This is illustrated in Figure 
8a, showing two closely adjacent vertices V x and 
being replaced by a single vertex V y The second pass 
checks again along the entire contour, but now for seg- 
ments with a length larger than the maximum length 
'max- ' n * n ' s case sucn an edge segment is divided into 
two shorter ones of equal length, as illustrated in Figure 
8b. In this Figure a vertex is inserted halfway in be- 
tween vertices Vj and necessarily the latter vertex 
is to be relabelled to V l+2 . The values of / min and / max are 
derived from the length / des . The relation between these 
parameters is constrained by not admitting an oscillato- 
ry behaviour, in which vertices are repeatedly removed 
in one resampling action and inserted again in the next. 
This leads to the requirement that / max > 2/ min . In the 
prototype implementation the / min was set at 1 /£/ des and 
/ max at 3/ min with satisfactory results. 
[0051 ] When a vertex is placed in the contour it is also 
assigned a velocity v, and acceleration a,- which are ob- 
tained by averaging the velocities and accelerations of 
the two vertices which are replaced (Figure 8a) or be- 
tween which it inserted (Figure 8b), which is necessary 
to maintain continuity in the dynamic situation. 
[0052] A prototype implementation of the described 
method was applied to clinical images, in particular to a 
typical clinical intravascular ultrasound image. The im- 
age represented a cross-section through a blood vessel, 
scanned by a full 2k radians rotation of a transducer 
mounted on the tip of a catheter. Such an image com- 
prises a number of more or less elliptical shapes, rep- 
resenting the layered structure of the vessel wall. Some 
of these elliptical shapes show up as bright ridges, while 
others show up as dark valleys. The area inside the ves- 
sel was partly blocked. The remaining open area was 
visible as a dark region. Both bright ridges and dark val- 
leys are interesting from a clinical point of view, as well 
as the edge of the open area. 

[0053] There are a number of reasons why it would 
be very difficult, if not impossible, to extract the desired 
contours from an image like this, using conventional (lo- 
cal) edge extraction methods. As usual in ultrasound im- 
ages, the signal to noise ratio is very low and large parts 
of the image are obscured by shadows. An additional 
difficulty is caused by the presence of the catheter (in 
the lower part of the open area, touching the vessel wall) 
and overlaid artificial tickmarks for distance measure- 
ment. The invented method succeeded well in defining 
the contours which describe the light and dark elliptical 
shapes in the vessel wall as well as the edge of the open 
area inside the vessel. 

[0054] For the bright ridges the inverse pixel grayval- 



ue was used as the image feature energy distribution, 
in order to allow the method to direct the contour to fol- 
low valleys. For the dark valleys the pixel grayvalues 
themselves can be used. 

[0055] In the contour to be defined along the edge of 
the open area inside the vessel, it was chosen to posi- 
tion the desired path of this contour between dark and 
bright image regions, so to make it to follow a path of 
maximum grayvalue gradient. As the contour in the used 
implementation of the method will follow valleys through 
the external image feature distribution, the reciprocal 
length of the grayvalue gradient is suitable as the image 
feature energy distribution. In such distribution, parts of 
the edge of the open area show up as valleys which do 
however not form a closed shape. The contour found by 
the method and the original image appeared to be an 
accurate definition of the edge of the open area. 
[0056] It appears that the operator-defined initial con- 
tour may have a very low resolution without jeopardizing 
the contour finding process. The initial resolution is au- 
tomatically and quickly increased by the resampling 
mechanism, which is integrated in the deformation proc- 
ess, until it is on the level that was specified via the pa- 
rameter / des . This will be realized after only a few defor- 
mation steps, and from there on the resampling process 
will keep the resolution at the required level. Conse- 
quently, the amount of user interaction in defining an in- 
itial contour can be very small. The amount of interaction 
may be reduced even further, by making use of the pos- 
sibility to create an open contour as an initial contour, 
and letting it grow by finding its own path. 
[0057] In case several contours have to be defined in 
more or less similar images, operator intervention can 
be reduced by using the resulting contour of an image 
as the start contour for a next image. This occurs, for 
example, in MRI or CT images in which a volume of a 
patient is scanned and the results are presented as a 
plurality of two-dimensional images, each image repre- 
senting a slice from the volume. 
[0058] In Figure 9 a workstation 50 is sketched in 
which the present invention can be applied. The work- 
station 50 comprises a display screen 51 for displaying 
images, a memory 52 and computing means 53. The 
workstation further comprises an input device 54 for in- 
put of data that represent an image in which contours 
are to be determined. Such an input device is, for ex- 
ample a magnetic or optical disc reader or a magnetic 
tape reader. It may also be a link to a central data stor- 
age or to a device that generates image data, for exam- 
ple an ultrasound device, a CT-scanner or an MRI-ap- 
paratus. The workstation, even can form an integral part 
of such device and allow the operator to determine con- 
tours in generated images and use the results obtained 
to determine which further image data should be taken. 
An operator console 55, for example a keyboard, is pro- 
vided for the user to give commands to the workstation. 
A pointing device 56, for example a lightpen, a tracker 
ball or a "mouse" for pointing at specific locations to de- 
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fine a start contour may be present as well. Additional 
devices may be connected to the workstation, such as 
a unit 57 to provide a hard copy of the image visible on 
the display or a storage unit 58 for storing images and 
contour data obtained for later recall. 
[0059] In the workstation data are read from the input 
device 54 and are stored in memory 52. The data rep- 
resent a parameter with a variation in a space in which 
contours are to be determined. The parameter data are 
normally given as a number of values, each value giving 
the parameter value as determined in a small element 
(pixel or voxel) of the space investigated. 



Claims 

1. Method for determining a contour (10) of an object 
in a space with a density variation, said density var- 
iation representing the object, comprising: 

defining a plurality of initial vertices (11,12,13) 
being located on a seed contour, said vertices 
being connected by edges (14, 15,16) to neigh- 
bouring vertices; 

starting from the seed contour, derive a dynam- 
ic contour that modifies its shape, 
defining an energy function comprising an in- 
ternal portion which is a function of the local 
curvature of the dynamic contour at each of the 
vertices and an external portion determined by 
the density variation at the location of each of 
said vertices; 

determining a final contour by continuously var- 
ying the position of the vertices till the energy 
function reaches a minimum, 

characterized in that for each of the vertices, 

the internal portion of the energy function is de- 
pendent on the angle between the edges con- 
necting the respective vertex to neighbouring 
vertices and in that 

the external portion of the energy function is re- 
lated to the density variation along a direction 
having equal angles with each of said edges 
connecting the respective vertex to neighbour- 
ing vertices. 

2. Method according to claim 1 , characterised in that 
the internal energy portion is derived from the dif- 
ference in angles of the edges at a series of adja- 
cent vertices. 

3. Method according to claim 2, characterised in that 
said internal energy portion is obtained by the con- 
volution of said difference in angles with a symmet- 
ric discrete filter having a zero frequency compo- 
nent equal to zero. 



4. Method according to claim 1 , 2 or 3, characterised 
in that the seed contour is open. 

5. Method according to claim 1 , 2, 3 or 4, character- 
5 ised in that if during minimization of the energy 

function an edge length exceeds a first predeter- 
mined threshold, the edge is replaced by an addi- 
tional vertex and connecting edges. 

10 6. Method according to claim 1, 2, 3, 4 or 5, charac- 
terised in that if during minimization of the energy 
function an edge length falls below a second pre- 
determined threshold, the edge and the connected 
vertices are replaced by a single vertex. 

15 

7. Method according to any one of the preceding 
claims, characterised in that a user defined ener- 
gy distribution is added to the external portion of the 
energy function. 

20 

8. Method according to claim 1 , characterised in that 
firstly a contour is determined in a first space and 
that, secondly, a contour is determined in a second 
space having a density variation having similar fea- 

25 tures to the density variation in the first space 
whereby the seed contour in the second space is 
defined by transfer of the final contour obtained in 
the first space. 



30 9. Arrangement (50) for determining a contour of an 
object in a space in which a given parameter has a 
value varying throughout said space, said parame- 
ter representing the object and the arrangement 
comprising 

35 

a memory (52) for storing parameter values in 
said space in the form of a matrix of pixel val- 
ues; 

a display (51 ) having a two-dimensional matrix 

40 of pixels; 

pointer means for indicating initial vertex posi- 
tions for the definition of a seed contour; 
computing means (53), having access to said 
memory (52) and said vertex positions, said 

45 computing means being adapted to 

derive a dynamic contour that modifies its 
shape starting from the seed contour,said com- 
puting means being adapted to 
evaluae an energy function comprising an in- 

50 ternal portion which is a function of the local 

curvature of the dynamic contour at each of the 
vertices and an external portion determined by 
the density variation at the location of each of 
said vertices; 

55 - said computing means being further adapted to 
determine a final contour by continuously var- 
ying the position of the vertices till the energy 
function reaches a minimum, said arrangement 
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being characterised in that 

for each of the vertices the internal portion of 
the energy function is dependent on the angle 
between the edges connecting the respective 
vertex to neighbouring vertices and the external 5 
portion of the energy function is related. to the 
density variation in a direction having equal an- 
gles with each of said edges connecting the re- 
spective vertex to neighbouring vertices. 



Patentanspruche 

1 . Verfahreh zur Bestimmung einer Kontur (1 0) eines 
Gegenstandes in einem Raum mit Dichteschwan- 15 
kungen, wobei die Dichteschwankungen den Ge- 
genstand darstellen, das folgendes umfasst: 



der interne Anteil der Energiefunktion von dem 
Winkel zwischen den den entsprechenden 
Scheitelpunkt mit den benachbarten Scheitel- 
punkten verbindenden Randern abhangt und 45 
dass 

der externe Anteil der Energiefunktion mit den 
Dichteschwankungen entlang einer Richtung 
verknupft ist, die den gleichen Winkel zu jedem 
der den entsprechenden Scheitelpunkt mit den 50 
benachbarten Scheitelpunkten verbindenden 
Rander hat. 

2. Verfahren nach Anspruch 1, dadurch gekenn- 
zeichnet, dass der interne Energieanteil von der 55 
Winkeldifferenz der Rander an einer Reihe benach- 
barter Scheitelpunkte abgeleitet wird. 



3. Verfahren nach Anspruch 2, dadurch gekenn- 
zeichnet, dass der interne Energieanteil durch Fal- 
tung der genannten Winkeldifferenz mittels eines 
symmetrischen, diskreten Filters mit einer Nullfre- 
quenzkomponente gieich Null erzielt wird. 

4. Verfahren nach den Anspruchen 1 , 2 oder 3, da- 
durch gekennzeichnet, dass die Kernelemente- 
Kontur often ist. 

5. Verfahren nach den Anspruchen 1 , 2, 3 Oder 4, da- 
durch gekennzeichnet, dass im Fall, dass wah- 
rend der Minimierung der Energiefunktion eine 
Randlange einen ersten vorher festgelegten 
Schwellenwert uberschreitet, der Rand durch einen 
zusatzlichen Scheitelpunkt und verbindende Ran- 
der ersetzt wird. 

6. Verfahren nach den Anspruchen 1, 2, 3, 4 oder 5, 
dadurch gekennzeichnet, dass im Fall, dass wah- 
rend der Minimierung der Energiefunktion eine 
Randlange einen zweiten vorher festgelegten 
Schwellenwert unterschreitet, der Rand und die mit 
ihm verbundenen Scheitelpunkte durch einen ein- 
zigen Scheitelpunkt ersetzt werden. 

7. Verfahren nach einem der vorhergehenden Anspru- 
che, dadurch gekennzeichnet, dass eine benut- 
zerdefinierte Energieverteilung zu dem externen 
Anteil der Energiefunktion hinzugefugt wird. 

8. Verfahren nach Anspruch 1, dadurch gekenn- 
zeichnet, dass erstens eine Kontur in einem ersten 
Raum bestimmt wird und zweitens eine Kontur in 
einem zweiten Raum bestimmt wird, die eine Dich- 
teschwankung mit den gleichen Merkmalen wie die 
Dichteschwankung im ersten Raum aufweist, wo- 
durch die Kernelemente-Kontur im zweiten Raum 
durch die Ubertragung der in dem ersten Raum er- 
haltenen, endgultigen Kontur definiert wird. 

9. Anordnung (50) zum Bestimmen einer Kontur eines 
Gegenstandes in einem Raum, in dem ein gegebe- 
ner Parameter einen Wert hat, der in dem genann- 
ten Raum variiert, wobei der Parameter den Gegen- 
stand darstellt und die Anordnung folgendes um- 
fasst: 

einen Speicher (52) zum Speichern von Para- 
meterwerten in dem genannten Raum in Form 
einer Matrix aus Pixelwerten; 
eine Anzeige (51 ) mit einer zweidimensionalen 
Matrix aus Pixeln; 

Zeigemittel zum Anzeigen von anfanglichen 
- Scheiteipunktpositionen zur Definition einer 
Kernelemente-Kontur; 

Rechenmittel (53), die Zugriff auf den genann- 
ten Speicher (52) und die genannten Scheitel- 



Definieren einer Vielzahl von Anfangsscheitel- 
punkten (11,12,13) auf einer Kernelemente- 20 
Kontur, wobei die genannten Scheitelpunkte 
durch Rander (14,15,16) mit benachbarten 
Scheitelpunkten verbunden sind; 
ausgehend von der Kernelemente-Kontur Ab- 
leiten einer dynamischen Kontur, die ihre Form 25 
verandert; 

Definieren einer Energiefunktion, die einen in- 
ternen Anteil, der eine Funktion der Krummung 
der dynamischen Kontur an jedem der Schei- 
telpunkte ist, und einen externen Anteil enthalt, 30 
der durch die Dichteschwankung an der Posi- 
tion jedes der genannten Scheitelpunkte be- 
stimmt wird; 

Bestimmen einer endgultigen Kontur durch 
standiges Andern der Position der Scheitel- 35 
punkte, bis die Energiefunktion ein Minimum 
erreicht, 

dadurch gekennzeichnet, dass fur jeden der 
Scheitelpunkte *o 
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punktpositionen haben, wobei die genannten 
Rechenmittel so ausgelegt sind, dass sie 
ausgehend von der Kernel emente-Kontur eine 
dynamische Konturableiten, die ihre Form ver- 
andert; 5 
eine Energiefunktion evaluieren, die einen in- 
ternen Anteil, der eine Funktion der lokalen 
Krummung der dynamischen Kontur an jedem 
der Scheitelpunkte ist, und einen externen An- 
teil enthalt, der durch die Dichteschwankung an 10 
der Position jedes der genannten Scheitel- 
punkte bestimmt wird; 

die genannten Rechenmittel ferner so ausge- 
legt sind, dass sie 

eine endgultige Kontur bestimmen, indem sie 15 
kontinuierlich die Position der Scheitelpunkte 
verandert, bis die Energiefunktion ein Minimum 
erreicht, wobei die genannte Anordnung da- 
durch gekennzeichnet ist, dass 
fur jeden der Scheitelpunkte der interne Anteil 20 
der Energiefunktion von dem Winkel zwischen 
den den entsprechenden Scheitelpunkt mit den 
benachbarten Scheitelpunkten verbindenden 
Randern abhangt und dass 
der externe Anteil der Energiefunktion mit den 25 
Dichteschwankungen entlang einer Richtung 
verkniipft ist, die den gleichen Winkel zu jedem 
der den entsprechenden Scheitelpunkt mit den 
benachbarten Scheitelpunkten verbindenden 
Rander hat. 30 



Revendications 

1 . Procede pour determiner un contour ( 1 0) d'un objet 35 
dans un espace presentant une variation de densi- 
te, ladite variation de densite representant I'objet, 
comprenant : 

la definition d'une plurality de sommets initiaux *o 
(11, 12, 13) etant situes sur un contour d'origi- 
ne, lesdits sommets etant relies par des bords 
(14, 15, 16) a des sommets voisins ; 
en partant du contour d'origine, la derivation 
d'un contour dynamique qui modifie sa forme ; 
la definition d'une fonction energie comprenant 
une partie interne qui est fonction de la courbu- 
re locale du contour dynamique au niveau de 
chacun des sommets et une partie externe de- 
terminee par la variation de densite a I'empla- so 
cement de chacun desdits sommets ; 
ia determination d'un contour final en variant de 
maniere continue la position des sommets jus- 
qu'a ce que la fonction energie atteigne un 
minimum ; 55 

caracterise en ce que, pour chacun des sommets, 



la partie interne de la fonction energie depend 
de Tangle entre les bords reliant le sommet res- 
pectif a des sommets voisins, et en ce que 
la partie externe de la fonction energie est Ii6e 
a la variation de densite suivant une direction 
presentant des angles egaux avec chacun des- 
dits bords reliant le sommet respectif a des 
sommets voisins. 

2. Procede suivant la revendication 1 , caracterise en 
ce que la partie d'energie interne est derived de la 
difference entre les angles des bords au niveau 
d'une serie de sommets adjacents. 

3. Procede suivant la revendication 2, caracterise en 
ce que ladite partie d'energie interne est obtenue 
par la convolution de ladite difference entre les an- 
gles avec un filtre discret symetrique presentant 
une composante frequencies nulle egale a zero. 

4. Procede suivant la revendication 1 , 2 ou 3, carac- 
terise en ce que le contour d'origine est ouvert. 

5. Procede suivant la revendication 1 , 2, 3 ou 4, ca- 
racterise en ce que si, pendant la minimisation de 
la fonction energie une longueur de bord d6passe 
un premier seuil predetermine, le bord est remplace 
par un sommet et des bords de liaison supplemen- 
taires. 

6. Procede suivant la revendication 1 , 2, 3, 4 ou 5, ca- 
racterise en ce que si, pendant la minimisation de 
la fonction energie, une longueur de bord tombe en 
dessous d'un deuxieme seuil predetermine, le bord 
et les sommets de liaison sont remplaces par un 
sommet unique. 

7. Procede suivant Tune quelconque des revendica- 
tions precedentes, caracterise en ce qu'une repar- 
tition d'energie definie par I'utilisateur est ajoutee a 
la partie externe de la fonction energie. 

8. Procede suivant la revendication 1 , caracterise en 
ce que, premierement, un contour est determine 
dans un premier espace et en ce que, deuxieme- 
ment, un contour est determine dans un deuxieme 
espace presentant une variation de densite com- 
portant des particularites semblables a la variation 
de densite dans le premier espace de telle sorte que 
le contour d'origine dans le deuxieme espace est 
defini par transfert du contour final obtenu dans le 
premier espace. 

9. Montage (50) pour determiner un contour d'un objet 
dans un espace dans lequel un parametre donne 
comporte une valeur variant sur Pensemble dudit 
espace, ledit parametre representant I'objet et le 
montage comprenant : 
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une memoire (52) pour stocker des vaieurs de 
parametre dans ledit espace sous la forme 
d'une matrice de vaieurs de pixel ; 
un affichage (51) comportant une matrice bidi- 
mensionnelle de pixels ; 5 
des moyens de pointeur pour indiquer des po- 
sitions de sommet initiales pour la definition 
d'un contour d'origine ; 

des moyens de calcul (53), ayant acces a ladite 
memoire (52) et auxdites positions de sommet, 10 
lesdits moyens de calcul etant propres a 
deriver un contour dynamique qui modifie sa 
forme en partant du contour d'origine, lesdits 
moyens de calcul etant propres a 
evaluer une fonction energie comprenant une 15 
partie interne qui est fonction de la courbure lo- 
cale du contour dynamique au niveau de cha- 
cun des sommets et une partie externe deter- 
minee par la variation de densite a I'emplace- 
ment de chacun desdits sommets ; 20 
lesdits moyens de calcul etant en outre propres 
a 

determiner un contour final en variant de ma- 
niere continue la position des sommets jusqu'a 
ce que la fonction energie atteigne un mini- 25 
mum, ledit montage etant caracterise en ce 
que 

pour chacun des sommets, la partie interne de 
la fonction energie depend de Tangle entre les 
bords reliant le sommet respectif a des som- 30 
mets voisins et la partie externe de la fonction 
energie est liee a la variation de densite dans 
une direction presentant des angles egaux 
avec chacun desdits bords reliant le sommet 
respectif a des sommets voisins. 35 
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FIG.1 
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